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1 Tethered Surfaces and Random Manifolds 

Several important developments of theoretical physics in the last 15 years 
come from the extension of the concept of random walk to fluctuating ex- 
tended objects. This has been very fruitful both in high energy physics, 
where the quantum fluctuations of strings (1 -|- 1-dimensional objects) and of 
p-branes {p + 1-dimensional) in Minkowski space are considered, and in con- 
densed matter physics, where the thermal fluctuations of 2-dimensional films 
or membranes in Euclidean 3-dimensional space are a fascinating subject (see 
for instance |]I[). 

It is known that a 2-dimensional surface, with an intrinsic metric Qa/s, 
embedded in flat c/-dimensional target space (the embedding being described 
by the mapping x f{x)) is characterized by its extrinsic metric 
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Dimensional analysis shows that the relevant terms involving r in the action 
may involve one derivative daT (tangent vectors) and two derivatives dapT 
(extrinsic curvatures) of the embedded surface. Two very different classes 
of models exist, characterized by the coupling between the intrinsic and the 
extrinsic metric. 



*Lectures given at the Cargese Summer School "Low Dimensional Applications of 
Quantum Field Theory", July 1995. 
t Physique Theorique CNRS 

■'■Laboratoire de la Direction des Sciences de la Matiere du CEA 



1 



(1) The intrinsic metric is proportional to the extrinsic one, i.e. is coupled to 
the embedding by the constraint gap = dardf^f. This is the case for the "rigid 
string" model of Polyakov, and for the equivalent Canham-Helfrich model of 
fluid membranes with bending rigidity. The action is 

S[r\ = Jd'xy^[T + ^{Af}'] (2) 

with r the bare string (or surface) tension, and k the bending rigidity. Renor- 
malization group calculations show that at large distances, the bending rigid- 
ity becomes irrelevant, and that the effective action is nothing but that of the 
Polyakov string, with an effective string tension t^s dynamically generated 

[|, i- 

(2) The intrinsic metric ga/s does not fluctuate (or it has a dynamics decoupled 
from that of the extrinsic metric). This class of models is not very useful 
for high energy physics, but is relevant in statistical physics to describe the 
so-called "tethered membranes", which generalize the concept of flexible 1- 
dimensional chains (polymers) to 2-dimensional networks It is this class 
of models that I shall discuss in these lectures. 

Let me consider a flexible two dimensional regular triangular network 
fluctuating in 3-dimensional space. A simple discrete action for the model is 

S = -K J2 nfUf + T Y {fi~ Tjf (3) 

neighbouring links 2, J 

triangles t,t' 

where fit is the normal vector to the triangle t. The first term is the ex- 
trinsic curvature term, with k the bending rigidity, the second term is a 
Gaussian elastic term. Numerical simulations and analytical arguments in- 
dicate a very interesting behavior for such a model. If k is small enough, 
the bending rigidity is irrelevant at large distances, and the surface is in a 
"crumpled", or "collapsed", phase, with {nt) = 0. This means that 0(3) 
rotational invariance is not broken and that the surface has no average ori- 
entation. If K is large enough, the surface is in a flat phase, characterized 
by a non-zero average orientation (n) ^ 0, and a spontaneous breakdown of 
0(3) invariance. These two phases are separated by a crumpling transition 
at some Kc- Numerical simulations indicate that (at d = 3) this transition is 
continuous, and characterized by non-trivial critical exponents. 
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An continuous model a la Landau- Ginzburg to describe tethered surfaces 
is given by the effective action 

S = Jd^x 

with hai3 = daTdpf the extrinsic metric. D is the internal dimension of 
the surface (the model describes in fact D-dimensional tethered manifolds), 
A = dada is the Laplacian, k the bending rigidity, t, K and /x are the effec- 
tive elastic moduli {t corresponds to a tension, K and are related to the 
so-called Lame coefficients). Neglecting the fluctuations, the minimization of 
this action shows that if t > 0, {daV) = and (ha/s) = 0, so that the surface 
is crumpled, while if t < (dar) 7^ and the surface is flat, (hajs) vanishes 
and there is a continuous crumpling transition at t = 0. This analysis ne- 
glects the effect of fluctuations and is valid only if the internal dimension 
of the manifold is D > 4. For D < A fluctuations become important, and 
their effect can be estimated by an e-expansion for D = 4 — e [^], or by a 
large d expansion P, ^ § (where d is the dimension of the target space). 
It might seems surprising that for 2-dimensional manifolds {D = 2) the flat 
phase phase still exists, since it is characterized by a spontaneous breakdown 
of the continuous rotational 0{d) symmetry, which should be forbidden by 
the Mermin- Wagner- Coleman theorem. Such a crumpling transition is in- 
deed forbidden for D = 1 (polymers): for any dimension > 1 of target 
space, infinite semi-flexible polymers with a non-zero bending rigidity k are 
always crumpled at large distance, and the flat phase does not exist. In fact 
there is no contradiction for D = 2. The transverse degrees of freedom of 
the manifolds r± (undulations) are coupled to the longitudinal degrees of 
freedom f\\ (phonons), so that the global symmetry of the model is not the 
compact group 0{d) (rotations), but the non-compact group of Euclidean 
displacements E((i) (translations+rotations). There is a non-trivial coupling 
between phonons (longitudinal modes) and undulations (transverse modes) 
which generates effective long range interactions between these transverse 
modes. In the presence of such long range interactions the Mermin- Wagner 
theorem does not apply. 

Let me give a tentative picture which emerge from the analytical and 
numerical studies of the crumpling transition. For D > 4 the crumpling 
transition is continuous and its critical exponents are given by mean field 
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Figure 1: Nature of the crumpling transition for phantom manifolds as a 
function of D and d 



theory. For D = 4 — e, e small, one loop calculations indicates that the 
crumpling transition is continuous for d > dc large enough, but becomes 
discontinuous (fluctuation induced first-order) for d < d^ small. Large d 
calculations shows that there is a second order crumpling transition for D > 
Di{d) < 2 with Di{oo) = 2. For D < Di{d) there is no fiat phase and 
no crumpling transition. Thus we expect that the domain 1 > D > 4 in 
the (c?, D) plane will be separated into three regions. A domain A for small 
D where the manifold is always crumpled, a domain B where there is a 
second order crumpling transition, and a domain C where the crumpling 
transition is always first order. Numerical simulations indicates that the 
point = 2, (i = 3) is in B Recent studies of the folding problem show 
that the point (D = 2, ci = 2) is in C 

The crumpled phase is the simplest to characterize. It corresponds to 
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t > in (^), and at large distance only the quadratic t term is relevant. 
Therefore the action for a crumpled manifold is Gaussian, and is nothing but 
the massless free field action 

S = Jd^'x^iV^f^' (5) 

The properties of such manifolds are easy to compute. For instance, if one 
considers a finite manifold with an internal extent L, its average squared size 
(_R^) in target space scales as 

r L(2-^) ifD<2 
oc <^ ln(L) ifD = 2 (6) 
constant ii D > 2 

which implies that the fractal dimension of a crumpled Gaussian manifold is 
df = 2D/{2 - D)iiD <2, and is infinite if D > 2. 

2 Self-avoiding crumpled Manifolds 

The above considerations apply to "phantom manifolds", which are free to 
intersect themselves. Indeed, the action (^) takes into account only local 
couplings in the internal space. Such local couplings are the only relevant 
one for strings, but for physical tethered networks self-avoiding interactions, 
which involve elements of the manifold which are close in target space but are 
arbitrarily far apart in internal space, are relevant. It is expected that such 
interactions will change the scaling properties of the manifold in the crumpled 
phase and at the crumpling transition. For instance, in the crumpled phase, 
the average squared size will now scale as 

o: L^" Q<u<l (7) 

with a critical exponent v > vq = sup[(2 — Z})/2,0]. Vi vq < v < 1 the 
manifold will be swollen by self-avoidance, but still crumpled. This is the 
case for polymers D = 1 for 1 < < 4, with for instance z/ = 3/4 for d = 2. 
The effect of self-avoidance may even be so strong that the manifold stays 
flat and that the crumpling transition disappear, in this case u = 1. Such a 
behavior has been observed in several numerical simulations of self-avoiding 
tethered membranes in three dimensions. 
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A simple analytical model to describe crumpled self-avoiding manifolds 



has been introduced in |TT|, |T2[. It is a simple extension of the continuous 
Edwards model for polymers. The action S (the free energy for a configu- 
ration r) is the sum of a Gaussian elastic energy and of a 2-body repulsive 
interaction, proportional to the coupling constant b: 

S[r\ = J d^x ^(V,f)2 + 6 Jd^xJ d^y 5'^ (rix) - f{y)) . (8) 

The internal dimension D may be taken as a continuous parameter, inter- 
polating between polymers {D = 1) and membranes {D = 2). The issue 
is to compute the critical exponents describing the scaling behavior of large 
manifolds, for instance the exponent u (related to the fractal dimension of 
the manifold by z/ = D/df), and the configuration exponent 7, related to the 
scaling of the partition function Z of a finite manifold with internal extent 
Lby 

Z (X L"' ^ constanf^ . (9) 

The mean field exponents are obtained by setting 6 = 0. One recovers the free 
Gaussian action and the exponents z/q = {2 — D)/2 and 70 = 1 — d{2 — D)/2 
(if D is not integer). 

Dimensional analysis shows that the mean field theory is invalid if the 
engineering dimension of b, e, is positive 

[b] = e = 2D- d{2 - D)/2 > . (10) 

In this case, we expect that u > uq. For small enough d, and certainly 
for d < D, we expect that the manifold is flat. The general picture of 
the expected behavior as a function of the internal dimension D and of the 
external dimension d is presented in Fig. ^ 

A natural idea is to compute the corrections to mean field by a e-expansion 
a la Wilson- Fisher. This has been done for polymers, first by using the so- 
called de Gennes trick [|1^: in the scaling limit the self-avoiding walk can 
be mapped (by a Legendre transform) onto a local field theory with 0(n) 
symmetry defined in the (i-dimensional target space, in the limit n — 0. 
It is then equivalent to study the limit of a single very long polymer and 
the massless limit of the n = theory, for which standard renormalization 
group theory is appliquable. Unfortunately no such equivalence exists for 
manifolds, beyond the case D = 1. Another renormalization scheme used 
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Figure 2: Self-avoiding manifolds as a function of d and D 



for polymers is the so-called direct renormalization scheme ||1 
applied to self-avoiding manifolds by Aronowitz & Lubensky 
Kardar & Nelson 
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It has been 
11| and by 

who performed calculations to first order in e. The 
basic idea of this method is to perform explicit perturbative calculations for 
a finite manifold and for e > 0. Perturbation theory is then UV and IR finite, 
but has UV divergences when e — 0. It appears that these poles in 1/e can be 
removed by reexpressing the observables in terms of adequate dimensionless 
renormalized quantities, such as the second virial coefficient. The internal 
size L of the manifold plays the role of the inverse of a renormalization mass 
scale, and renormalization group equations can be obtained by considering 
the L dependence of the renormalized theory. At first order in perturbation 
theory, the consistency of these renormalization group equations has been 
checked explicitly by Duplantier, Hwa & Kardar |T^. However, at that 
time it was not clear whether this direct renormalization approach could 
be justified beyond first order (except for D = 1, where the de Gennes 
trick is used to show the equivalence between direct renormalization and the 
standard minimal subtraction scheme). 
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3 Renormalization for multi-local Theories: 



Recently is became possible to prove the consistency of this approach, and 
the renormalizability of the model (H) directly in the internal D-dimensional 
space, despites the fact that this model is a non-local field theory in D 
dimensions [0. Let me present the general idea for the proof, which is due 
to B. Duplantier, E. Guitter and myself. The perturbation theory for this 
model is obtained by expanding the observables as power series in h. The 
bi-local "interaction vertex" (in field theoretic language I call it a bilocal 
operator) is written in Fourier transform as 

, , = 6\r{x) - r{y)) = [ d% e^^(^»-^^(^)) . (11) 



It can be viewed in a Coulomb gas representation as the integral over the 
"charge" A; of a neutral "dipole" with charge +k at x and charge —k at y. The 
term of order in the perturbative expansion of the partition function (as 
well as of other observables) involves K dipoles (xi,?/i), . . . , {xk.Vk)- The 
integration over the charges fci, . . . kK gives an integral over the positions of 
the dipoles of the determinant of the "dipole energy" quadratic form Q 

j ■■■ jWd^x^ d% det [Q[xi, Vi]] ""^ . (12) 
1=1 

K ^ 

Q is a K X K matrix such that J2 kiQijkj is the Coulomb energy (in D- 
dimensions) of the K dipoles. Each Qij is a linear combination 

Qij = Go{xi, Xj) + Go{yi, yj) - Go{xi, yj) - Go{xj, yi) . (13) 

of the Coulomb potentials Go between the endpoints of pairs of dipoles i and 
3 

Go{x,x') = {rix)rix% = ^^^^^j^ \x - xf-^ • (14) 

This Coulomb potential is properly defined for < D < 2 by analytic con- 
tinuation in D. For < D < 2 it is negative, but it vanishes for x = x' (while 
for Z) > 2 it diverges), and it decreases at large distances (as for D > 2). 
The integration over the 2K points in a non-integer D-dimensional space 
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Figure 3: A non- neutral regular configuration and two neutral UV singular 
configurations 

can also be defined properly by analytic continuation in D and the use of 
distance geometry. This amounts to replace the integration over the 2K x D 
coordinates of the 2K points by an integration over the K x {2K — 1) scalar 
distances between these points. 

One can show that when the determinant det[(5] vanishes short distance 
UV singularities occur in the integrals. This occurs if and only if some of the 
end-points of (not necessarily the same) dipoles coincide, so that the end- 
points form "atoms" , while the dipoles form "molecules" , and if moreover one 
can assign non-zero charges ki to the dipoles while each atom stays globally 
neutral. This condition is more easily depicted graphically on Figure 3. 

The associated singularities of these integrals are related to the behavior 
at short distance of the expectation value (with respects to the free Gaussian 



model) of products of bifocal operators as given by Equ. 11. One can show 
that this short distance behavior is encoded in a multilocal operator product 
expansion (MOPE) , which generalizes Wilson's operator product expansion. 
Let me give two examples: 

When the two points x and y of the bi-local interaction operator tend 
towards a single point, this operator can be expanded in terms of local oper- 
ators involving derivatives of the field r. The first terms of the expansion are 
explicitly (not writing explicitly the D and d dependence of the coefficients) 




^~xy +ci\x- y\'-^-^{x'' - i/")(x^ - /) : Vaf V/jf : 

+■■■ (15) 

1 is the identity operator (its expectation value is 1), the : : in the operator 
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: V^r V/3r : denotes the normal ordering subtraction prescription required to 
deal properly with the UV singularities contained in V^r V/jr. 

The second example is less simple, and shows that when the end-points 
of two bilocal operators tend pairwise towards two different points, this gen- 
erates again bilocal operators 



^^'T '-^'J' = do 



-d/2 



^ ^' ^'^2 + ... (16) 

This structure is generic, and products of local and bilocal operators 
generate multilocal operators of the general form 

<^{x,,---,xp} = fd% n [i^^ror S''{ro-r{x^)) A{x,)] . (17) 
•' i=i 

where the Ai{xi) are local operators, which can be decomposed into products 
of multiple x-derivatives of r. The rrii are integers. For P = 1 and m = 
one recovers local operators A{x) (m > gives 0). For P = 2, mi = 1712 = 
and Ai = A2 = 1 one recovers the bilocal interaction operator, etc. . . These 
operators have a very special form: they can be viewed as a local convolution 
in the target d-dimensional r space of a non-local product (in the internal 
D-dimensional space) of the P local operators Ai. 

The MOPE implies that the formalism of renormalization theory and of 
renormalization group equations, which has been developed for local quan- 
tum field theories, can be adapted for this model. One is in fact interested 
in the IR scaling behavior of the lattice model, when some length scale L 
goes to 00. This lattice model is described by the action (|), with a short 
distance lattice cut-off a. To study this IR limit it is equivalent to look at the 
UV continuum limit of the model when the physical length scale L is kept 
fixed, while the UV cut-off a goes to 0. In this limit one can construct, via 
renormalization, a finite renormalized theory with a = 0, which obeys renor- 
malization group equations. From these equations, one recovers the large 
distance behavior of the lattice model we started from. The procedure works 
well in perturbation theory when one is close to the upper critical dimension, 
i.e. for e small, and it leads to the e-expansion. 

In our case, the MOPE can be used to determine, by power counting, 
which multilocal operators are relevant and give UV singularities (poles in 
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Figure 5: Examples of nested singular configurations 



1/e). Then one can also show that these poles can be subtracted by adding 
to the action counterterms proportional to the marginally relevant multi- 
local operators, leading to the UV finite renormalized theory. For the model 
of self-avoiding manifolds, this analysis shows that the UV divergences are 
associated only with local and bilocal operators, as depicted on Fig. ^ and 
that only three operators are relevant: the identity operator 1, the elastic 
energy operator (Vr)^ and the bilocal operator 6'^{f{x) —f{y)). 1 is strongly 
relevant, and gives power-like UV divergences proportional to {a being 
a short-distance cut-off). The two other operators are superficially relevant, 
they give logarithmic UV divergences or equivalently poles in 1/e at e = 0. 

The fact that the so-called superficial divergences, associated to a global 
contraction of points towards a singular configuration, can be subtracted by 
counterterms is a consequence of the MOPE. A complete proof of the renor- 
malizability of the theory is possible, but much more delicate. It requires a 
control of the sub divergences coming from successive contractions associated 
to nested singular configurations, such as those depicted of Fig. ^ . 



11 



Figure 6: The /5-f unction and the RG flow for e > 0. 



4 Scaling for infinite self-avoiding Manifold 

A first apphcation of this formalism is the derivation of scaling laws. Since 
the model is renormalizable (at least perturbatively), it can be made UV 
finite (for e ~ 0) by introducing two counterterms in the action. The new 
renormalized action is of the form 

5[r1 = ^ Jd^'x (Vf)2 + bnfi' Z, J J d'^xd^ 5{r{x) - r{y)) . (18) 

6ij is the dimensionless renormalized coupling constant (the perturbative ex- 
pansion in hji is UV finite order by order). Z is a wave-function renormal- 
ization factor and Z}, a coupling constant renormalization factor, both are 
perturbative series in with poles up to degree l/e^~^ at order K. /i is 
the renormalization momentum scale. As for ordinary local theories, such as 
the Landau- Ginzburg- Wilson $^ model, one can change hR and r in Equ. |TB| 
into bare quantities in order to rewrite the renormalized Hamiltonian as a 
bare action given by Equ. |. The renormalization group /^-function and the 
anomalous dimension 7 of the field r are defined in the standard way 



d 



7(M = -If^^^-z 



(19) 
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The /3-function is found to be of the form 

Pibn) = -ebji + cbji + 0{bj^) ; c = c(D) positive constant , (20) 

and therefore there is, at least for small e > 0, an IR attractive fixed point 
6* = 0{e), which governs the scaling behavior of self-avoiding polymerized 
surfaces at large distance. The existence of this fixed point ensures the uni- 
versality of this non-trivial scaling for e > 0, and that no new interactions, 
possibly non-local in external space, are generated by the RG transforma- 
tions. 

The explicit calculation for the scaling exponents u and 7 leads to the 
same results for the scaling exponents u and 7 than the direct renormalization 
method at first order in e. With this method higher order calculations are 
feasible, but technically quite difficult. In particular, already at second order 
the RG functions cannot be expressed analytically, and numerical integration 
methods have to be developed. Work is in progress to compute the scaling 
exponents at order e^. 

5 Finite Size Scaling and direct Renormaliza- 
tion 

The model given by Equ. § describes an infinite manifold with fiat internal 
metric, corresponding to an infinite and regular fiexible lattice. Finite man- 
ifolds are described by a similar model, but the D-dimensional manifold M 
is now embodied with a fixed non-trivial Riemannian metric gap{x) (exam- 
ples are closed manifolds with the topology of the sphere iSd or the torus 
T'^), and may have a boundary dM (open manifold with the topology of 
the disk for instance). A similar analysis can be performed for such models, 
and the MOPE structure of short distance singularities is still valid, but new 
local operators A(x), which depend on the internal metric on M and on the 
boundary dM, such as the scalar curvature i?, appear in the MOPE and in 
Equ. (0). The renormalized action now contains at least five operators and 
five independent renormalization factors Z 

S[r\ = f Zl + f Z{Vrf + f f Zb6'^{f-r) 

+ [ ZR + f Zl . (21) 

Jm JdM 
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The curvature operator R is superficially relevant only for D = 2 and the 
boundary operator Jgj^ 1 only for D = 1. When these additional terms are 
not relevant, the first three renormalization factors Z are the same for finite 
curved manifolds than for the infinite fiat plane. This property is analogous to 
the renormalization property of local field theories in finite geometries, which 
justifies the finite scaling laws for critical systems in finite geometries, and it 
has two very important consequences: (i) The scaling hypothesis at the basis 
of the direct renormalization approach, which relies explicitly on calculations 
with finite manifolds, is shown to be valid to all orders in perturbation theory; 
(ii) for "abstract" manifolds with dimension D < 2, with the only exception 
of open polymers {D = 1 open surface), the following hyperscaling relation 
[O relating the configuration and the u exponents holds: 

= 1 - ud . (22) 



6 Self-avoiding Manifold at the tricritical 0- 
point: 

Finally, let me briefly discuss recent results obtained with K. Wiese on the 
scaling behavior of polymerized membranes at the 0-point [jl8|. This point 



separates the swollen phase, where the self-avoidance repulsive forces that I 
considered previously dominate, from the dense collapsed phase, where short 
ranged attractive forces dominate. At the 9-point the effective two body 
repulsive coupling b vanishes, and two different interactions may become 
relevant. The first one is the 3-body contact repulsion, which is usually 
considered for polymers 



c/^x rfi/ d'^z S^rix) - r{y)) 5%f{x) - f{z)) . (23) 



The second one is a modified 2-body interaction, repulsive at short range but 
attractive at larger range (A^ is the ci-dimensional Laplacian) 



{i • =- I I d^xd^y Ar- 5" {f{x)-f{y)) . (24) 
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Figure 7: Relevant interactions for the 0-point in the d-D plane. 



Calculations at first order are not feasible analytically, and already require 
numerical evaluations of complicated integrals. The results of such one loop 
calculations are schematically depicted on Fig. ^, where the domains where 
the 3-body and modified 2-body terms are respectively relevant are shown. 
This indicates that the last modified 2-body term is the relevant one for 
2-dimensional manifold in any external dimension d. There is also a quite 
interesting and non-trivial crossover between the two terms around D = 4/3 
(i = 6, which must be studied by a double e-expansion. 



7 Conclusion: 

The theoretical study of the scaling behavior of polymerized flexible mem- 
branes leads to the development of new multilocal continuum field theories, 
and to new applications of renormalization group methods. I hope that these 
methods will lead to a quantitative progress in the understanding of the be- 
havior of real 2-dimensional polymerized membranes. This requires results 
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beyond first order in the e-expansion (recall that D = 2 correspond to e = 4), 
and a better understanding of the relation between this RG approach and 
more heuristic or approximate methods, such as variational methods or ap- 
proximate recursion relations. The sophisticated renormalization theory for 
multilocal models presented here should hopefully also find applications in 
other problems of statistical physics, or in other areas of theoretical physics. 
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